Segmental relaxation in semicrystalline polymers: a mean field model for the 
distribution of relaxation times in confined regimes 
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The effect of confinement in the segmental relaxation of polymers is considered. On the basis of a 
thermodynamic model we discuss the emerging relevance of the fast degrees of freedom in stimulating 
the much slower segmental relaxation, as an effect of the constraints at the walls of the amorphous 
regions. In the case that confinement is due to the presence of crystalline domains, a quasi-poissonian 
distribution of local constraining conditions is derived as a result of thermodynamic equilibrium. 
This implies that the average free energy barrier AF for conformational rearrangement is of the 
same order of the dispersion of the barrier heights, S(AF), around AF. As an example, we apply 
the results to the analysis of the a-relaxation as observed by dielectric broad band spectroscopy in 
semicrystalline poly(ethylene terephthalate) cold-crystallized from either an isotropic or an oriented 
glass. It is found that in the latter case the regions of cooperative rearrangement are significantly 
larger than in the former. 



INTRODUCTION 

Segmental relaxation is a process whereby a certain 
initial chain conformation, produced by a thermody- 
namic fluctuation in an amorphous region of a polymer, 
evolves towards a more probable state. This mechanism 
has a central role in several processes, such as elastic 
recovery 0L dielectric response diffusion of small 
molecules [3], glass transition 0] and local chain read- 
justments at the boundaries of growing crystalline do- 
mains [5|, to mention a few. 

Recently, the problem of understanding how con- 
finement affects conformational relaxation has attracted 
growing interest [f|. In this respect, semicrystalline poly- 
mers represent a natural ground where some progress 
may be done in this field, since neither melt- nor cold- 
crystallization lead to complete crystalline states in these 
materials. At the end of the crystallization process, the 
remaining amorphous regions cannot convert to crys- 
talline because of the conformational constraints induced 
by the confining crystal domains through pinning and 
chain connectivity. In particular, the dynamics of the a- 
relaxation associated to the segmental motions above the 
glass transition temperature T g is strongly affected by the 
progressive development of the crystalline phase @, S, 0] • 

From the thermodynamic point of view, the establish- 
ment of constraints in the amorphous regions reduces the 
conformational entropy, causing an increase of the chem- 
ical potential /i a with respect to the value fi a o of the 
unconstrained liquid. On the other hand, crystallization 
from an undercooled liquid (either constrained or not) is 
associated to a decrease of the chemical potential down 
to a value /i c < fi a o- 

The crystals reduce chain mobility; this is not only 
due to the decrease of the volume accessible to the amor- 
phous chains, but also because many of these chains are 
anchored to the crystals. Thus, crystal growth is always 



accompanied by an increase of fi a ■ In the initial stages of 
the crystallization process, the difference fj-a—l^ao is small 
compared to the chemical potential drop Afi ac = \i c — /i a o 
associated to the amorphous-to-crystalline transition of a 
monomer. However, as crystallinity increases, the change 
of (i a — /i a Q accompanying further transition of monomers 
to the crystalline state is so large that it cannot be com- 
pensated by the corresponding free energy decrease as- 
sociated to \Afi ac \. The whole system thus tends to a 
stationary state beyond which the phase transition can- 
not proceed further. 

On the basis of this idea, the vitrification process ob- 
served in the low-temperature crystallization of uniax- 
ially oriented poly(ethylene terephthalate) (PET) has 
been recently discussed (lfj| . 

Experiments show that the glass transition tempera- 
ture T g increases with respect to a completely amorphous 
material, as an effect of a finite crystallinity [HI, H3 |. 
However, it is important to point out that this varia- 
tion is always accompanied by an increase of the temper- 
ature interval over which the amorphous component of 
the system de- vitrifies upon heating. These effects, which 
can be easily observed by means of differential scanning 
calorimetry, are particularly evident when a dynamical 
mechanical thermal analysis in the a-relaxation region 
is performed (see e.g. ref. jlO] for the case of PET). 
This suggests that the emergence of constraints to chain 
readjustments and the development of inhomogeneities 
in their character could be inherently connected. 

The analysis of the crystallization process by means 
of broad band dielectric spectroscopy (BDS) gives sim- 
ilar indications. Recent investigations on the crystal- 
lization of poly(propylene succinate) and poly(propylene 
adipate) (PPS and PPA respectively) [bj, [lj] represent 
particularly clear examples of the broadening effect as- 
sociated to the increase of constrainings. These stud- 
ies pointed out that the dominant characteristics of the 
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a relaxation change during crystallization. Namely, the 
rather fast and comparatively narrow dielectric a-process 
observable in completely amorphous samples is progres- 
sively substituted by a much slower and broader process 
(in frequency) eventually overwhelming the former. In 
some polymers like PPA and PET [3, 0, EH , the two 
relaxations can be observed as distinct, separable pro- 
cesses, one increasing in intensity at the expenses of the 
other as crystallization proceeds. In other polymers, like 
PPS 0, Ej and poly(L lactide acid) [H El El one 
observes just one process whose average characteristic 
relaxation time decreases progressively while broadening 
in frequency. 

With the aim of describing how the establishment of 
structural constraints hinders conformational rearrange- 
ment and gives rise to time heterogeneity, we shall resort 
to the concept of cooperatively rearranging region (CRR) 
originally introduced by Adam and Gibbs for the descrip- 
tion of liquids close to their glass transition [2(| . 

A CRR is the smallest amorphous domain where a con- 
formational rearrangement may occur without causing 
any structural change at its boundary. Last condition is 
in fact a constraint which determines the existence of a 
free energy barrier AF that a thermodynamic fluctua- 
tion needs to overcome to produce a conformational re- 
arrangement EH0]- This, on turn, affects the charac- 
teristic relaxation time of the CRR through 



,AF/k B T 



(1) 



where ks is the Boltzmann constant and T is the absolute 
temperature. 

In a similar manner we shall consider small amorphous 
regions embedded by crystals. In general, however, AF 
does not only depend on the conditions at the bound- 
ary of the confining volume, but also on the actual chain 
conformation within the CRR. In other words, once a 
conformational transition has occurred, the free energy 
barrier in the new state in general differs from its former 
value. This already suggests a possible origin of the con- 
nection between the mean barrier AF and its dispersion 
5(AF) in an ensemble of CRR's. 

The description of the evolution dynamics of AF in a 
CRR is a difficult task. For this reason we shall limit 
ourselves to propose a statistical mechanical model with 
the aim of pointing out the characteristic features that an 
ensemble of CRR's in a semicrystalline polymer should 
have in order to meet the general principles of thermody- 
namics. This will lead us to the derivation of a distribu- 
tion of constraining conditions, i.e., of single time decay 
processes contributing to the overall segmental relaxation 
process. 

As an example of application, this distribution will be 
used to analyze some dielectric relaxation data of PET 
already available in the literature. 
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FIG. 1: Schematic of an amorphous domain viewed as a 
thermodynamic subsystem. In the presence of confinement by 
crystals, the direct interaction of the conformational degrees 
of freedom with the heat bath vanishes. 



PHENOMENOLOGICAL MODEL 
Thermodynamic considerations 

A finite amorphous domain in a polymer is a thermo- 
dynamic subsystem which, from the statistical mechani- 
cal point of view, can be thought to be made of two parts; 
one consists of only the chain conformational degrees of 
freedom (C-component), while the other (V-component) 
collects the vibrational ones. The observation that the 
specific heat of a glass is very close to that of the cor- 
responding crystal (wherever the system is able to crys- 
tallize) suggests that these two components are approxi- 
mately independent. In other words, 



b Z CO nf - 



(2) 



where Z, Z V ib and Z con f are the partition functions of 
the whole subsystem, the vibrational component and the 
conformational one respectively. However, there are pro- 
cesses where the direct C-V interaction is important and 
approximations like eq. [2] must be handled with care. 
Segmental relaxation in confined regimes is one of these 
cases. 

The scheme reported in Fig. 1 represents our subsys- 
tem, with its two components C and V, interacting with 
the heat bath (HB). The V component of the subsys- 
tem always interacts directly with both the HB and the 
conformational part C. On the other hand, the direct 
coupling of the latter with the HB is strongly affected by 
confinement. When the subsystem is embedded by crys- 
tals, for instance, the conformational degrees of freedom 
are isolated from the HB, the conformational fluctuation 
modes with wavelength larger than the characteristic di- 
mensions of the subsystem are depressed and the C-HB 
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coupling is rendered possible only through the C-V inter- 
action. Hence, the latter is responsible of the fluctuations 
occurring in the conformational component under crystal 
confinement. 

Of course the above one is an idealized picture. In the 
case of polymers, vibrational and conformational degrees 
of freedom are not so sharply separated with regards to 
confinement effects. Rather, we have to think that there 
are motions that are sensitive to boundary conditions, 
others that are not, and a relatively wide region of the 
spectrum in between. We can safely say that the modes 
for which confinement is uneffective are comparatively 
faster than the others, but we cannot exclude that suit- 
ably fast segmental modes belong to the intermediate 
class. This depends of course on how effective is the con- 
finement and on the degree of cooperativity characteris- 
tic of the system. We shall comment this point further 
below. 

As illustrated above, the effect of the crystals in regime 
of confined relaxation is the quenching of the conforma- 
tional modes which are associated to density fluctuations 
with wavelengths larger than a certain value. In this re- 
spect, it is worth mentioning that there exists another 
situation where the long wavelenght fluctuation modes 
appear depressed, namely, when their characteristic evo- 
lution time is so slow to appear frozen with respect to 
the observation time scale (this is in fact the mechanism 
involved in the glass transition). When the attempt is 
made of describing the apparent relaxation dynamics in 
terms of statistical thermodynamics, we are necessarily 
lead to consider the existence of virtual conformational 
constraints. The difference with respect to confinement 
by crystal consists in the absence of a localized volume, 
with ideal rigid walls, which could cause the quenching 
of the modes. In other words, the constraints are spread 
out over the entire liquid, as they are only of a kinetic 
nature. For this reason we call this mechanism kinetic 
confinement as opposed to the former which we shall call 
static. 

In the case of kinetic confinement the direct C-HB cou- 
pling is depressed by the quenching of some interaction 
channels (i.e. the long wavelength modes), but does not 
totally disappear. In this situation the direct C-HB and 
V-C interactions could be of comparable relevance in af- 
fecting the relaxation dynamics. 

In the following we shall consider cooperatively rear- 
ranging regions whose confinement is mainly of a static 
nature, in the sense that wherever the border of these 
regions is not coincident with the surface of a crystal, 
the heat transfer through this "hole" is negligible with 
respect to the transfer via the V-C coupling. For con- 
venience we shall call any of them static-CRR (SCRR). 
The case of the kinetic confinement regime is left to fu- 
ture work. 



The mean field distribution of constraints 

Consider a (small) domain of volume V in an uncon- 
strained amorphous region of a polymer, and let Z u be 
the partition function associated to it. The emergence of 
conformational constraints in this domain has two main 
consequences. On the one hand, the partition function 
decreases to Z' u < Z u because some states are ruled 
out; on the other, the transition from one conformational 
state to another in general involves the passage through 
non-equilibrium states, which means that a fluctuation 
free energy barrier AF > has to be overcome in or- 
der that the conformational change takes place. From 
the statistical mechanical point of view, the set of these 
non-equilibrium states is described by a reduced parti- 
tion function, Z c , collecting only the states (among all 
those contributing to Z' u ) which allow for conformational 
rearrangements [20(: 



AF 



-k B T\n-^-. 



(3) 



The deviation from an initial (metastable) equilibrium 
state of a SCRR - a subsystem in contact with a heat 
bath - is in general accompanied by an increase of its 
energy (see e.g. ref. [21|, Sect. 21). Within the spirit of 
a mean-field description of the conformational transition 
process, we follow the suggestion indicated by ref. [iol ] 
and take Z c as the part of Z' u consisting of just the terms 
with energy E > E m i n , being E m i n a threshold value 
below which conformational rearrangements are not pos- 
sible. Note that the establishment of an E m i n is related 
to the transition process only. Of course, there may be 
cases where the actual E m i n and the actual energy value 
Ei of the initial conformational state coincide, but on av- 
erage the relation E m i n > Ei holds. In our mean field 
treatment the energy threshold E m i n will be taken as a 
measure of the level of constraints. 

We now consider an ensemble of SCRR's with a certain 
distribution of shapes and dimensions, that is, the num- 
ber of monomers within the SCRR's is distributed around 
an effective mean z. For convenience, in the following, we 
shall refer the extensive quantities to the monomers - as 
the smallest crystallizable units - and deal with chemical 
potentials instead of free energies. 

In connection with the cooperativity character of con- 
formational rearrangement, all monomers in a SCRR 
are equivalent, so we can associate to each of them 
a specific partition function Z = [Zscrr] 1 ^ , where 
Zscrr = Zscrr(z) is the partition function associ- 
ated to the whole SCRR. The specific partition function 
Z collects the complexity of the monomer viewed as a 
thermodynamic subsystem; furthermore, it implicitly in- 
cludes the interaction with the other monomers forming 
the SCRR. Associating to each monomer a specific par- 
tition function, has the same meaning as associating to 
it a chemical potential (cf. ref. 20]). 
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In general, when a conformational change occurs in a 
SCRR, the new state differs from the parent one also 
in the value of the lower energy bound per monomer, 
C, separating the fluctuations that can induce rearrange- 
ment from the others which cannot. With regards to the 
(instantaneous) characteristic relaxation time tscrr, the 
energy ( suitably defines the state of a SCRR or, which is 
the same, of all the monomers of which is made; the rea- 
son is that different chain conformations yield the same 
t SCRR provided the associated £ values are the same. 

We now distribute all monomers into different classes, 
labelled by the value of C characterizing the SCRR's to 
which they belong. At this point the individuality of the 
SCRR's is lost and we remain with the problem of finding 
the most probable monomers' distribution p((). 

It is clear that any change in the distribution is ren- 
dered possible by the monomers which undergo a fluctu- 
ation to a state with an energy larger than the threshold 
C labelling the class to which they belong; let's call "ac- 
tive" these monomers. Once this threshold is exceeded, 
a monomer may change its state and fall into a different 
class; this is to say that the SCRR to which it belongs 
has reached a different conformational state. We cannot 
say which class an active monomer will fall in after each 
transition (it can even remain in the same class), but in a 
stationary state the ensemble of all active monomers will 
be distributed among the different classes in the most 
probable way. 

In order to find this distribution, we consider the whole 
set of active monomers as a thermodynamic subsystem 
in thermal equilibrium with the heat bath. The total 
volume of this subsystem is constant (it is related to the 
overall crystallinity) , so the most probable distribution 
can be found by minimizing the associated Helmholtz 
potential 



reduced specific partition function 



A = U - TS . 



(4) 



where S is the entropy and U is the average energy in- 
volved in the transfer of monomers from one class to the 
other. More explicitly, we associate to each class the av- 
erage energy {E}^ of the corresponding active monomers; 
then, U is a suitable average of these terms (see below). 

On one hand, the entropy term takes the customary 
form for a distribution of units in the different classes 



S = -k B J2p(0 lnp(C) 

C 



(5) 



on the other, the average energy involved in the transfer 
of units between classes is 



= $>(CK(0 (£> c 
c 



(6) 



where the average energy (E)^ , associated to the 
monomers able to leave their class (£), is related to the 



-E m /k B T 



m; E m >Q 

(being E m the energy of the m-th state) by 
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the factor 



{E) ^~ d(ksT)-^ nZ ^ ] 



Z 



(7) 



(8) 



(9) 



with Zq = (( = 0) , is the probability that a fluctuation 
increases the monomer energy above the actual minimum 
threshold for a rearrangement. 

The equilibrium distribution minimizes the Helmholtz 
potential eq. |4] under the condition 



= ^2p(C) A/i(C) = const. , 

C 

where the chemical potential barrier is given by 



A/M s ; ' ( |j 



(10) 



(11) 



This condition simply says that, independent of how the 
exchange of elements among classes occurs, the thermo- 
dynamic state of the liquid, i.e. the average level of 
conformational constraints (to which the rearrangement 
chemical potential barrier height A/i is associated) is sta- 
tionary, at least within the time scale of the observation. 
Shouldn't this condition being fulfilled, we would envis- 
age the possibile variation of the crystallinity during the 
observation. 

Straightforward calculations yield the following form 
for the most probable distribution: 



p(0 = iV(A) e -^(?)^c/ fc - T - AA ^) 



(12) 



where N(X) is the normalization constant and A, the La- 
grange multiplier associated to the condition eq. [101 is 
related to the average chemical potential barrier A/z by 
the implicit relation 



Afj, 



3 1 



-w r (0(E)Jk B T-XAfi 



C>0 



(13) 



The fact that the probability factor w r is not contained 
in the entropy term may appear inconsistent; however, 
this is not the case. The point is that any monomer what- 
soever, independent of whether it will change its state or 
not, is an element of a class; for this reason the entropy 
must be written in the above form eq. El On the other 
hand, not all the elements of a class undergo a change 
of state to another class, but only the fraction w r (0- 
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To see this more clearly, consider that the equilibrium 
distribution could be as well obtained by extremizing S 
with the condition on Ap,, eq. [TU1 together with the con- 
ditions of constancy (stationarity) for the average energy 
involved in the exchange of elements among the classes, 
i.e. U — const. The expression of the average energy, 
of course, has to take into account that not all the ele- 
ments of a class [oc p(C)] undergo a transition, but only 
the fraction w r (£) of them. 

It is easy to show that also in the framework of this 
thermodynamics of constraints, the variance of A/x can 
be expressed in terms of the mean value Ap by means of 
the expression 



-a— 2 5A M 



[S(Ap)f ee <A M 2 ) - An = 



OX 



(14) 



This equation is important, because it is related with the 
observations reported in the introductory section con- 
cerning the connection between free energy barrier height 
and heterogeneity, both introduced by the constraints. 
When the A An term dominates in the exponent of eq. [T2l 
the distribution approaches a poissonian form. This is to 
say that the square root of the variance and the mean 
An value tend to be the same in this limit. 



Relaxation function 

Segmental motion underlies the fluctuation regres- 
sion of quantities like, for example, density or elec- 
tric polarization. A key quantity for the description 
of these processes is the relaxation function <j>(t) = 
(B(t)B(O)) /[B(0)] 2 , that is, the normalized autocorre- 
lation function of a time dependent quantity B{t) with 
respect to its value at the initial time t = 0. 

In general, cooperative rearrangement characterizes 
the segmental relaxation dynamics at long times [22j |. 
However, in the case of static confinement, where the 
kinetic component in the conformational constraints is 
small, the cooperative character of the motion is expected 
to extend down to comparatively shorter times. Consid- 
ering the confined segmental motion only, we represent 
the relaxation function as a superposition of single-time 
decay contributions: 



t 



-(C) 



where the characteristic relaxation time, 
r(0 



= T * e AF(()/k B T ^ 



(15) 



(16) 



carries information on cooperativity through the free en- 
ergy barrier 



associated to the rearrangements of the SCRR's, and r* 
is a phenomenological triggering time which is expected 
to be rather smaller than the central relaxation time To 
of the a-process. Indeed, we showed at the beginning 
of this section that, if static confinement is effective, the 
conformational transitions must be stimulated by much 
faster (local) process. 

The above expression of <fi(t) still remains undeter- 
mined unless is known. On the other hand, the cal- 
culation of this partition function is in general not trivial 
at all and requires further modeling (which is out of the 
scope of the present contribution). 

Alternatively, with the aim to illustrate a practical ap- 
plication of the model, we express the partition function 
in the form 



/oc 
dEp(E,n)e- E/kBT ee Z c , 



(18) 



where the density of states p(E, n) accounts for degener- 
acy and amplitude of the energy intervals between suc- 
cessive levels and it is taken in the form 



p(E,n)=p E n , 



(19) 



with n an integer to be determined by fitting, and po an 
irrelevant constant (cf. eqs.l8ll9land[TTj). In fact, n selects 
the relevant energy interval of the states contributing to 
Z' u . Note that the above density of states, eq. HH is the 
same as for a set of n+l independent harmonic oscillators 
in contact with a heat bath [23|. 

In the case that n = the free energy barrier AF takes 
the form 



AF(C) ^z(. 



(20) 



AF(C) = zA/i(C) 



(17) 



We recall that this functional dependence has been 
adopted to treat the case of strong confinement (i.e. the 
present static confinement) in ref. 24j on the basis of 
qualitative considerations. Equation [20] states that the 
free energy of cooperative rearrangement has no entropy 
contribution. While this may be reasonable at high ( 
values, it is not an acceptable approximation when £ is 
small, i.e., when the unconstrained liquid state is ap- 
proached. 



EXAMPLE OF APPLICATION 

In this section we present the analysis of dielectric 
relaxation data as an example of application of the 
above model. We shall consider the cases of PET cold- 
crystallized from either a nematic or an isotropic glass. 
The data have been already published and the details of 
the experimental procedures, together with some aspects 
of the data analysis, can be retrieved from the litera- 
ture [H, 24 1 . We only recall that, as a first step, the 
complex dielectric permittivity e — e' — ie" is measured 
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FIG. 2: Imaginary part of the dielectric permittivity at the 
temperature T m£a a = 96 °C, resolved in components, of an 
unoriented PET sample cold-crystallized at T c = 96 °C for 
360 min. The squares are the experimental values 



FIG. 3: Imaginary part of the dielectric permittivity at the 
temperature T meas = 130 °C, resolved in components, of a 
cold-drawn PET sample cold-crystallized at T c = 140 °C for 
120 min. The squares are the experimental values. 



as a function of the angular frequency lu. Then, the data 
are fitted to an empirical expression where each relax- 
ation process is described by a corresponding Havriliak- 
Negami term, i.e. 



A £j 



[1 + (iujTjp] 



(21) 



in which j labels the relaxation process, Asj is the 
strength of the relaxation, Tj is the central relaxation 
time, cij S ] 0, 1] describes the width of the relaxation 
(a 3 = 1 for the narrowest) and bj S ] 0, 1] its asymme- 
try (bj = 1 for a symmetric relaxation). The relaxation 
function associated to the j-th process is then obtained 
from the corresponding Havriliak-Negami component by 
means of a Fourier cosine transform }25j . i.e. 



/(*) = - 

7T 



cos(u)t) — 



I) 



Ae 7 - 



(22) 



In the following examples, we shall consider the a- 
process emerging in the presence of a finite crystallinity. 
This relaxation has been named a' in ref. [161 ] (it cor- 
responds to the a s process in ref. Q) and is believed 
to be a manifestation of static confinement. The faster 
a-relaxation is active also in the absence of crystallinity 
and, wherever present in the crystallized system, is rather 
associated to kinetic confinement and will not be dis- 
cussed here. 

Figure 2 reports the imaginary part of the dielectric 
permittivity e", as a function of the frequency v = u)/2ir, 
of an initially isotropic PET sample cold-crystallized for 
360 min at a temperature T c = 96 °C [l6j |. The different 
components derived from the fitting procedure are also 



shown, namely, the two distinct a processes, the tail of 
the /3-relaxation and the conductivity contribution. 

With respect to the a relaxation, the a' component is 
broader and centred at a lower frequency; this indicates 
that the latter process is influenced by constraining con- 
ditions to a larger extent than the former. 

Figure 3 shows the s" vs. v dependence at a tempera- 
ture T — 130 °C in the case of a cold-drawn PET sample 
cold-crystallized at T c — 140 °C for 120 min. Comparing 
it with the relaxation spectrum of the initially isotropic 
sample, we observe that the a process is absent and that 
the a' contribution is even broader than the correspond- 
ing one in Fig. 2. The central relaxation times of the two 
a' components reported in Figs. 2 and 3 are similar, that 
is, at the same T the a' process in oriented PET is much 
slower than in initially isotropic PET. 

Figure 4 reports in a log-log plot the relaxation func- 
tions (open squares and triangles) associated to the a' 
processes of Figs. 2 and 3, obtained after eq. fJU (the in- 
set shows the same in a lin-log plot). The solid lines are 
the best fitting results obtained by using the function 
defined by eq. 1151 which in a more explicit form reads: 



<f> a ,{t)=N{\n) d( 



exp 



Z, 



C,n+1 



Zo,n 
Z(,n 

T* \Zo,n 



\k B T 



4> a >(\,T*,z,n;t) 



(23) 



where the fitting parameters have been explicitly indi- 
cated in last row, Zq^ is defined by eq. [18] above and 
the term corresponding to the average energy in the ex- 
ponent of eq. [12] has been expressed as u>£,n {E) ^ n — 
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FIG. 4: Relaxation functions (f>(t) associated to the static 
constrained a-processes of PET reported in Figs. 2 and 3, as 
obtained from eg. 1221 The solid lines are the best fittings from 
cq. [23] The inset shows the same using linear ordinates. 



Z(,n+i/Zo in ; note that the normalization factor N now 
depends on n (cf. eq. [T2|) . 

The analysis is carried out by first setting a time lower 
limit, say t m im and fitting only the data for t > t m i n ; n 
is assigned also a starting value, say n = 0, which is set 
fixed during each fitting run. Then, the other parameters 
(A, r* and z) are let free to adjust. The procedure then 
consists in lowering t m i n progressively (while trying dif- 
ferent integers n) so to reach the minimum possible value 
compatible with an acceptable x 2 (of course, the condi- 
tion t* < t min must always hold). The reliability of the 
results can be assessed by checking if the approximate 
equality (cf. eq. QI 



fc^Tln 



zA/i 



(24) 



is verified by the final values of the fitting parameters. 

We report in Table 1 the values of crystallization tem- 
perature T c , measurement temperature T meas , width pa- 
rameter a (eq. l2"Tj) of the relaxation (in all cases b = 1) 
and central relaxation time To, together with the fitting 
parameters A, r*, z and n worked out by fitting c/)(t) 
for t > t m i n ; in the same table are also reported t m i ni 
the average energy threshold £, the dispersion 5( around 
it, the average chemical potential A/i and the dispersion 
around it, <5[A/i] (last quantities being calculated as or- 
dinary averages using the distribution p(£) once A and n 
have been found). 

The values reported in the Table verify eq. [M] above. 
Note that a more effective cooperativity (larger z) leads 
to smaller r* values at fixed ro and A/i (see Table I). 

There are a few points that deserve consideration. 
First of all, we always find that t* «C tq; this is con- 



sistent with the idea that the observed a' relaxation is 
indeed stimulated by faster processes. Furthermore, we 
can see that in the case of oriented PET the number of 
monomers of a SCRR is significantly larger than in the 
case of initially isotropic PET, although £ has a lower 
value (cf. Table I). This is an important result, because 
it is consistent with the observation that the vitrification 
of amorphous domains induced by (cold-) crystallization 
is indeed rather more pronounced in oriented than in 
isotropic PET [10]. Cold-drawing was indeed shown to 
affect the density autocorrelation length of PET, £, in 
the sense of increasing its value along the orientation 
direction, i.e. £„ em ~ 2^ so [1(3]. The present results 
thus suggest a very clear connection between cooperativ- 
ity and £. The fact that £ is a decreasing function of the 
temperature 1^, 26} further supports this results, since 

Zi so |t=96°C < 2„em|T=130°C- 

About the actual z values worked out by fitting, two 
aspects have to be taken into account. The first is that 
they come out from a mean field model, which means 
that they are not expected to be exact. The second is 
that they represent an underestimate when a (weak) di- 
rect coupling between the SCRR conformational degrees 
of freedom and the heat bath is effective: faster fluctu- 
ation decays are attributed to smaller SCRR's, within 
the model, because relaxation mechanisms other than 
through the local degrees of freedom are ruled out. 



We conclude this section by comparing the present 
results with those of our previous analysis on oriented 
241 ] . In that article we guessed that a gaussian 



PET 

distribution of constraining conditions could be a rea- 
sonable description of time heterogeneity; moreover, we 
approximated the barrier free energy AF like in eq. 1201 
above. The resulting fitting function was rather simpler 
than that used here for the analysis of the data. How- 
ever, the direct derivation of z and AF from the fitting 
procedure was not possible, although the latter quan- 
tity could be estimated from r* and to assuming that a 
relationship like eq. [24] would hold. We thus found that 
ksTh^TQ/T*) ~ 11.9 kcal/mol [3], to be compared with 
the present value of 8.9 kcal/mol. 



On the other hand, the estimate of the free energy 
barrier dispersion from the fitting of the relaxation func- 
tion, as reported in ref. [13], yielded 5(AF) ~ 8.8 
kcal/mol, while in the present case we obtain z5[A/j] ~ 
6.7 kcal/mol. This is to say that the former approach, 
ref. [2^], can be safely used for an estimate of 5(AF) in 
the static confinement regime; more noticeably, already 
that simple analysis indicated that S(AF) and AF were 
comparable, a circumstance that has been verified to oc- 
cur also in other semicrystalline polymers [27j |. 
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CONCLUDING REMARKS 

A central result of this paper is the connection between 
the average free energy barrier AF = zAfi associated 
to the conformational rearrangement within the SCRR's 
and the width SF = z5(Afi) of the barrier height distri- 
bution in an ensemble of SCRR's: 



AF ~ S(AF) . 



(25) 



In the introductory section, it was already pointed out 
that the increase of the barrier height in conformational 
rearrangement is accompanied by an increase in the het- 
erogeneity of constraining conditions. Furthermore, from 



an analysis of the relaxation behavior of PET in ref. [24 1 
it was found that a relationship like eq. [55] could hold. 
The theoretical model developed above represents an at- 
tempt to interpret this phenomenology. 

The classical theory of thermodynamic fluctuations (cf. 
e.g. ref. [2l[) already establishes a connection between 
the average value of an extensive quantity and the am- 
plitude of its fluctuations: as the volume of a subsystem 
decreases, the latter approaches the former. Analogously, 
the application of the present model to the analysis of 
the relaxation function, indicates that a relationship like 
eq. [^characterizes the fluctuations of constraints on sub- 
systems consisting of z monomers, which in fact turn out 
to be small. 

In the present model, the heterogeneity of the SCRR 
ensemble arises as a consequence of thermodynamic equi- 
librium wherever static confinement is effective. Of 
course, the question whether conformational relaxation 
under confinement (cither static or dynamic) consists of 
a superposition of single time decay processes, or it is the 
manifestation of an intrinsically non exponential decay, is 
not answered presently. This is just because of the ther- 
modynamic nature of the model: a framework where the 
dynamics of monomer transfers among different (^-classes 
is disregarded. 

We now turn to an aspect that we pointed out al- 
ready in the thermodynamics subsection, namely, the 
separation between relaxation modes that are sensitive 
or non-sensitive to the confining conditions. Comparing 
the analyses carried out above on the two PET samples, 
we observe that, although the central relaxation times 
To are practically the same (i.e. 5 and 4.2 xlO -3 sec), 
the fitting triggering times r* are significantly different; 
indeed, oriented PET is characterized by an enhanced co- 
operativity and a larger readjustment free energy barrier. 
On the other hand, the meaning of eqs.[T5landriBlis that 
the a! relaxation is the long time scale manifestation of 
fast processes through the coupling of the local degrees of 
freedom with the conformational ones. In a sense, there 
is no primary process defining for instance the central 
relaxation time To, because the latter just characterizes 
an asymptotic dynamics. Through the calculation of the 



free energy AF, our analysis ultimately indicates how 
faster than To can be a process for it being still affected 
by confinement (cf. eq. |24|) . Below this lower bound (t*) 
there may well be, depending on the cooperativity of the 
system, fast segmental modes which are local enough not 
to be affected by the presence of the crystals. 

A number of open points still remain to go into thor- 
oughly after this preliminar description of the restricted 
conformational dynamics, but probably the most impor- 
tant one would be the extension of the model to the 
case of the ordinary a-process, where confinement has 
a kinetic character. We hope, however, that the present 
contribution may be a useful starting point for further 
progress in this field. 
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Sample 


T c 


Tmeas 


a 


TO 


train 


A 


r* z 


n 


C 8( An 6[Afj] 




(°C) 


(°C) 




(sec) 


(sec) 


(kcal/mol)- 1 


(sec) 




(kcal/mol) (kcal/mol) (kcal/mol) (kcal/mol) 


isotropic 


96 


96 


0.4 


5 x 10"" 


3 x 10" 5 


1.45 


2 x 10" 5 2.5 


5 


6.5 2.1 1.65 1 


nematic 


140 


130 


0.18 


4.2 x 10" 3 


10" 6 


1.92 


5.4 x 10" 8 7.6 


5 


5.7 2.3 1.23 0.88 



TABLE I: Crystallization temperature T c , measurement tem- 
perature Tmeas, width parameter a, central relaxation time 
to (cf. refs. fwl , HJ), lower fitting time bound t m i n and 
relaxation function parameters (A, r*, z and n) of the ot- 
processes in static confinement regime for PET samples crys- 
tallized from isotropic and nematic glasses; the average min- 
imum energy threshold its dispersion around the mean SC,, 
the average chemical potential A/i and its dispersion <5[A/i] as 
calculated from the fitting parameters using the distribution 
eq. [12] are also reported. 



